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ABSTRACT 

Context. The “central engine” of AGN is thought to be powered by accretion on a central nucleus believed to be a super-massive black hole. 
The localization and exact mechanism of the energy release in AGN are still not well understood. 

Aims. We present observational evidence for the link between variability of the radio emission of the compact jet, optical and X-ray continua 
emission and ejections of new jet components in the radio galaxy 3C 390.3. 

Methods. The time delays between the light curves of the individual jet components and the light curve of the optical continuum are estimated 
by using minimization methods and the discret correlation function. 

Results. We find that the variations of the optical continuum are correlated with radio emission from a stationary feature in the jet. This 
correlation indicates that the source of variable non-thermal continuum radiation is located in the innermost part of the relativistic jet. 
Conclusions. We suggest that the continuum emission from the jet and counterjet ionizes material in a subrelativistic outflow surrounding the 
jet, which results in a formation of two conical regions with broad emission lines (in addition to the conventional broad line region around the 
central nucleus) at a distance S 0.4 parsecs from the central engine. Implications for modeling of the broad-line regions are discussed. 
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1. Introduction 


The variable continuum flux in AGN, signaling the activ¬ 
ity of the central engine, is detected throughout the entire 
electromagnetic spect r um, on time-s c ales from days to years 
JPeterson^t_^ l2002t IZhend Il996l: IWamsteker etZI Il997t 
IShapovalova et all 2001 ). The continuum flux is believed to 
be responsible for ionizing the cloud material in the broad¬ 
line region (BLR). Localization of the source of the variable 
continuum emission in AGN is therefore instrumental for un¬ 
derstanding the mechanism for release and transport of en¬ 
ergy in active galaxies. In radio-quiet AGN, representing about 
90 % of the AGN population, the presence of rapid X-ray flux 
variations and iron emission line (Fe Ka) indicates that most 
of the soft X-ray emiss ion originates from the accretion disk 
( iMushotzkv et alJl993h . In radio-loud AGN, the activity of the 
central engine is accompanied by highly-relativistic collimated 
outflows (jets) of plasma mat erial formed and accelerated in 
the vicinity of the black hole llFerrar lE^. Inhomogeneities 
in the jet plasma appear as a series of compact radio knots 
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(jet components) observed o n scales ranging from several light 
week s to about a kiloparsec jAlef et al.ll996l:lKellermann et alJ 
Continuum emission from the relativistic jet d ominates 
at all energies (see lUlrich et alJll997t IWorrallir2005l) . swamp¬ 
ing the X-ray emission associated with the accretion flow. 
Hence, the continuum variability in radio-loud AGN may be 
related to both the jet and the instabilitie s of accretion flows 
jMushotzkv et alJl993Uuirich et alJl997l) near the c entral en¬ 
gine. The unification scheme Jllrrv & Padovanill99.5h of radio- 
loud AGN suggests that the central powerful optical continuum 
and broad emission lines are viewed directly in radio quasars 
and BL Lacs, whereas in radio galaxies these can be hidden by 
an obscuring, dusty torus, and therefore, the bulk of continuum 
and broad-line emission in radio galaxies may be attributed to 
the relativistic jet rather than the central engine. The presence 
of a positive correlation between beamed synchrotron emis¬ 
sion from the base of the jet and optical nuclear emission in 
radio galaxies suggests that the optical e mission is non-thermal 
and may originate from a relativist ic jet Jchiaberge et al ■Il999l 
l2002t iHardcastle & Worrallll200(il) . The detection of a corre¬ 
lation between the radio and optical emission variability from 
the nuclear region would be the most direct evidence of opti- 
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Fig. 1. Radio structure of 3C 390.3 observed in 1998.21 with 
very long baseline interferometry at 15 GHz (2 cm). Innermost 
fraction of the jet is shown in the inset. The resolving point- 
spread function (beam) plotted in the upper left corner is 
0.87 X 0.55 mas oriented at an angle of 8.0 degrees (clockwise 
rotation). The peak flux density in the image is 190mJy/beam 
(3.1 X lO'^K) and the rms noise is 0.2mJy/beam (IJy = 
10 26Wm^^Hz '). The contours are drawn at 1, V2, 2... of 
the lowest contour shown at 0.6 mJy/beam. The structure ob¬ 
served is quantified by a set of two-dimensional, circular 
Gaussian features (shaded circ les) obtained fr om fitting the vis¬ 
ibility amplitudes and phases jPearso nll9?7l) . Similar fits have 
been obtaine d for ten observations fro m the 15 GHz VLBA sur¬ 
vey database ( iKellermann et alj2004t) . for the purpose of cross- 
identifying and tracing different features in the jet. The labels 
mark three stationary features (D, SI, and S2) and a subset of 
moving components (C2-C6) identified in the jet. Note that 
two more components, C7 and C8, have been first identified in 
the jet in the VLBA images at later epochs (see Fig.|2j. 

cal continuum emission coming from the jet. No observational 
evidence has yet been reported for a link between optical/UV 
continuum variability and the radio jet in radio-loud AGN. In 
Section 2, we analyse the structure and kinematics of the pc- 
scale jet in 3C 390.3 and look for correlations between variable 
emission of jet components and nuclear optical emission on 
scales less than one parsec. The structure and emission mecha¬ 
nism of the sub-parsec-scale region around the central nucleus 
is discussed in Section 3. In Section 4, we discuss the results 
and draw conclusions. 

2. The link between radio jet and optical 
continuum in the radio galaxy 3C 390.3 

To search for a relation between variability of the contin¬ 
uum flux and changes in the radio structure in AGN on 
scales of ~ Ipc (Ipc = 3.086 x 10^®m), we combine the 
results from monitoring of the radio-loud broad emission- 



Fig. 2. Separation of the jet components relative to the station¬ 
ary feature D (filled diamonds) for ten epochs of VLBI ob¬ 
servations. Moving components are denoted by leftward tri¬ 
angles (C4), rightward triangles (C5), upward triangles (C6), 
downward triangles (C7), leftward triangles (C8), and station¬ 
ary components SI and S2 are marked by filled circles. The 
lines represent the best linear least-squares fits to the compo¬ 
nent separations. 

line galaxy 3C 390.3 (re d shift z - 0.0561) in the optical 
ovalova et al.[ Eo^ ISergeev et al 1 l2nn2h . UV (|a end 
, and X-ray dLeighlv et al.lh997h regimes with ten very 
long baseline interfe rometry (VLBI) observ ations of its radio 
emission at 15 GHz dKellermann et al.ll2004l) made from 1992 
to 2002 using the VLBA'. 

2.L Structure and kinematics of the parsec-scaie jet 

To parameterize the structure of th e radio emissi on, we ap¬ 
plied the technique of modelfitting dPearso ^1 19971) and fit in¬ 
terferometric visibilities from each of the ten VLBA datasets 
by a set of circular Gaussian components (Fig. ^ and Tabled 
in the Appendix) and used the positions and flux densities of 
these components for tracing the evolution of the jet emis¬ 
sion on angular scales of ~ 3 milliarcseconds (mas). At the 
distance of 3C 390.3, 1 mas corresponds to a linear distance 
of 1.09 pc for the flat cosmology with the Hubble constant 
Hq - 70km s * Mpc * and the matter density Qn, = 0.3. 

From the component separations measured within 3 mas 
from the feature D at the narrow end of the jet (Fig. 0, we 
identified five moving components C4-C8 (in addition to the 
previously known components C2 and C3; Alef et al. 1996) 
and two stationary features SI and S2 separated from D by 
0.28 + 0.03 mas and 1.50 + 0.12 mas, respectively (Fig.|2j. The 
more distant stationary feature S2 can be explained by a small 
change in the direction of the flow that causes the relativistic 
dimming of the radio emission ilGomez et al.lll99^ . The pres¬ 
ence of SI may be related to both geometrical and physical 
factors. Linear fits to the observed separations of the compo- 

’ Very Long Baseline Array of National Radio Astronomy 
Observatory, Socorro, NM, USA 
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Fig. 3. The vari ations of the continuum fluxes of 3C 390.3 during the 1992 to 2002.2 time period. Optical con tinuum light 
curve at SIOOA dshanovalova et al]l200lU.Sergeev et alJl2002li (black squares), UV continuum fluxes ( fZhendfToOhh scaled by a 
factor of 50 (grey squares), and variations of radio flux density from D and S1 components (filled red diamonds and circles) are 
presented. The times of ejection, fo, of radio knots from D (presumed base of the jet) and the times of their separation, fsi, from 
the component S1 are marked by open triangles and filled triangles respectively. 


nent C4-C8 from the component D yield proper motions of 
0.2-0.4 mas/yr, which correspond to apparent speeds of 0.8- 
1.5 c, where c is the speed of light. Back-extrapolation of the 
fits to C5 and C6 indicates that these two components may 
originate from a single event. They could result from a mov¬ 
ing perturb ation in the jet that creates a forward and a reverse 
shock pair dOomez et al.lll99^ . We used the linear fits to es¬ 
timate, for each moving component, the epoch, to, at which it 
was ejected from the component D and the epoch, fsi, when 
it passed through the location of the stationary feature S1 (see 
Figs. 121 and|3i. 

2.2. Correlations between variable emission of jet 
components and optical continuum 

In Fig. 12] the properties of the components D and S1 iden¬ 
tified in the radio jet of 3C 390.3 are compared with the 
optical continuum emis sion at 5100 A l| Sergeev et al.l 1200^ 
Shapovalova et '^1200ih and 1850A llZhenJl996h . The radio 
fluxes from D (fo) and SI (/si) features (Fig.|3} show an ap¬ 


parent anticorrelation during 1996.4-2002.2 time period, which 
indicates that there should be a time delay between radio light 
curves. We should expect that the radio flare propagates down¬ 
stream of the jet, and hence the D component should be the 
leading source of emission. 

A poor radio sampling (10 points) makes it impossible to 
use the discrete correlation or cross-correlation functions for 
determining the time delay td-si between fo and /si (hereafter 
in the subscripts of r the first term indicates the leading source). 
Here, we employ the minimization method for linearly in¬ 
terpolated radio light curves. Each light curve is linearly inter¬ 
polated, and the data points of one light curve are shifted. For 
every value of the shift we calculate the normalized;^^^/(A-n), 
where N is the number of overlapping points and n - 2 is the 
number of fitted parameters. Application of this method shows 
that the best time delay is Td-si ~ 1 yr (Fig.|2} for both the en¬ 
tire radio data and for a subsample of eight radio observations 
made during the time period of optical measurements (from 
1994.6 to 2000.6). The time delay of ~ 1 yr is comparable with 
the mean time fo-si = 1-04 + 0.16 yr at which the moving 
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Fig. 4. The normalized statistics as a function of time delay 
between radio light-curves (ten data points between 1994.6 and 
2002.2) of D and S1 components. This dependence is presented 
for the entire radio data (dashed line) and for a subsample of 
eight radio epochs between 1994.6 and 2000 (solid line). The 
minimum corresponds to the best-fit delay. The best-fit delays 
are found in the range between 0.9 years to 1 year. 



Fig. 5. The significance level of the Spearman’s rank correla¬ 
tion coefficient for different time delays between /d and /opt 
(dashed line) and /opt and /d (solid line) light curves. The hor¬ 
izontal dotted lines correspond to the 99 % and 99.9 % confi¬ 
dence levels. 

components C4-C7 (the data for C8 are too sparse for making 
reliable estimates) pass the distance between D and SI. This 
indicates that the time delay between variations of radio light 
curves is related to the passage of moving components between 
D and SI stationary components. 

We use a different method to check for correlations between 
variability of the jet radio emission and the optical continuum 
emission at 5100A (/opt). Using the advantage of a dense sam¬ 
pling of optical measurements we generate paired radio-optical 
data points. For every value of time shift and for every radio 
measurement we assign an optical flux which is linearly in¬ 
terpolated between two nearest epochs of optical observations 
around the shifted epoch of a radio measurement. The signif¬ 



Fig. 6. The significance level of the Spearman’s rank correla¬ 
tion coefficient for different time delays between /si and /opt 
light curves. The horizontal dotted lines correspond to the 95 % 
and 99 % confidence levels. 


icance level of the Spearman’s rank correlation coefficient is 
calculated for different values of the shift. 

This method has been applied to the optical (106 points) 
and radio (eight points) data measured between 1992 and 2002. 
Correlations have been searched between the optical contin¬ 
uum and radio emission from the components D, SI, S2, and 
from a composite flux from D-i-Sl. No significant correlation 
exists for the moving features, the stationary feature S2 and the 
flux density of D-i-Sl (not shown in Fig.|3l. A significant cor¬ 
relation (> 99 %, see Fig. |3 is present between /d and /opt for 
time lags ranging from ro-opt ~ (1-2 - 1.6) yr. No significant 
correlation is found in the /opt - /d relation plane where the op¬ 
tical continuum is assumed to lead the radio emission (Fig.|3- 
This result strongly suggests that the radio emission precedes 
the optical continuum by approximately 1.4 yr, and that there 
should be a physical relation between variable radio emission 
from the feature D of the jet and optical continuum emission. 
One should expect that the optical continuum also correlates 
with and follows the radio emission of S1 component with time 
delay of Tsi-opt = ro-opt - td-si ~ 0.4 yr. Fig. shows that 
there is a positive correlation (at the S 95 % confidence level) 
in the /si - /opt relation plane for a wide range of time lags 
from ~ 0 yr to ~ 0.7 yr. To check this result we use the discrete 
correlation function jEdelson & Kroliljll998l) which yields a 
strong positive correlation between the radio flux density of 
SI (leading source) and the optical continuum for time lags 
ranging between 0 yr and 0.8 yr with slightly larger correlation 
function at small time delays, ~ (0 - 0.4) yr (see Fig.^}. Both 
statistical tests show a positive correlation between variations 
of /si and /opt for rsi-opt ~ (0 - 0.7) yr. It is most probable that 
the optical continuum is produced near or at the location of the 
radio emission (0-0.4 light year from the feature SI). If the op¬ 
tical continuum trails the radio emission then the optical and 
high-energy photons are likely to be produced by Compton up- 
scattering of radio p hotons by the relativistic electrons of the jet 
jUlrich et al.lll99'% . In either case, the size of the region from 
which the optical continuum is emitted should be very small. 
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Fig. 7. The discrete correlation function for different lags be¬ 
tween /si and /opt during the 1994.5 to 2000 time period. The 
0.1 year bin width of the time delay is adopted to minimize the 
errors of correlation function. Icr error bars are presented. 

This can also be deduced from the facts that, in 3C 390.3, the 
changes of the optical flux at 5100 A follow the variations in the 
UV a nd soft X-ray emissi on, with time delays of about 5 + 5 
days dPietrich et al.lll998i) and the overall spectrum in the op¬ 
tical to soft X-ray range is rep resented by a single pow er-law 
with a spectral index a - 0.89 dWamsteker et alJll997l) . These 
observations indicate that a large fraction of the variable con¬ 
tinuum emission is non-thermal and it is produced in a region 
of ~ 5 light days (0.004 pc) in size. 

The link between the optical continuum and the component 
SI in the radio jet is further supported by the correlation be¬ 
tween the continuum light curve and the characteristic epochs 
fo and fsi of the moving components C4-C7. For the compo¬ 
nents C4-C7, the epochs fsi of separation from the stationary 
feature SI are coincident, within the errors, with the maxima 
in the optical continuum (Fig.|3- All four ejection events oc¬ 
cur within ~ 0.3 yr after a local maximum is reached in the 
the intensity of the optical continuum (the average time delay 
between the maxima and the epochs fsi is 0.18 + 0.06 years). 
The null hypothesis that this happens by chance is rejected at 
a confidence level of 99.98%. This suggests that radio ejection 
events of the jet components are coupled with the long-term 
variability of optical continuum. 

3. The central sub-parsec-scale region 

The link between the radio emission from D and S1 suggests 
that the appearance of these stationary features in the jet is due 
to physical, rather than geometrical, factors. The location of 
stationary features with respect to the central nucleus is im¬ 
portant for understanding the structure of the central engine 
and its radiation mechanism producing the variable continuum 
emission. If the component D is identified with the proverbial 
“base” of the counterjet then the component S1 should be the 
base of the jet. We should expect that /d < fsi because of rela¬ 
tivistic Doppler effect. The fact that /d ^ /si during the six year 
time period (Fig.|3 rules out the assumption that D is the base 



Fig. 8. The variatio ns of soft X-ray fluxes at 0.1-2 keV 
iLeighlv et alJ ll 19971) for 3C 390.3 between 1995 and 1996 su¬ 
perimposed with the epochs of origin of C5 and C6. Icr error 
bars are presented for all data. 

of the counterjet. The component D is likely to be the base of 
the VLBI jet near the central engine; the base being the legation 
at which the flow either becomes supersonic ()Dalv et alJl988li 
or optically thin (iKoni 7illl98lHTohanovll99ar or releases the 
energ y contained in the Povnting flux jRomanova & Lovelac^ 
119961) . The identification of the component D with the cen¬ 
tral engine is supported by the correlation between the ejec¬ 
tion epoch, fo, of the component C5 and the variability of the 
X-ray flux at 0.1-2 keV (Fig.|8}. The ejection of C5 occurred 
after a dip in the X-r ay emission and hardening of the spectrum 
dTeighlv et al.lll997l) suggesting (similar to 3C 120; Marscher 
et al. 2002) that the soft X-ray emitting disk material disap¬ 
pears into the black hole and a fraction of the infalling matter 
is ejected into the jet. The component D should then be associ¬ 
ated with the radio emission coming from the immediate vicin¬ 
ity of the central engine of 3C 390.3 (assuming that absorption 
along the line of sight passing through outer layers of the torus 
is sufficiently weak). This emission is probably generated in or 
abov e a hot corona at a distance > 2 00 above the accretion 
disk (lFabianll2004t IPonti et alJl2004l) (Rg - 2GMhhfc^ is the 
Schwarzschild radius for a black hole of mass Mbh, where G is 
the Newtonian gravitational constant). 

The component SI then can be associated with the sta¬ 
tionary radio feature which may be produced by the inter- 
nal oblique shock f ormed in the continuous relativistic flow 
dOomez et alJll99^ . In this scenario, the energy release in the 
form of variable radiation can be produced by interaction of 
relativistically movin g compact compone nt or shock (?) with 
the stationary shocks dOomez. et al .111 99.51) produced at the po¬ 
sitions of D and SI of the jet. The time delay between variations 
of /d and /si is equal to the time (td-si = fo-si) at which the 
ejected component passes from D to S1. This is in fact the case 
for the 3C 390.3 (1 yr=s 1.04 yr, see the previous section). The 
characteristics of a long-term variability of optical continuum 
flux (amplitude of optical flares and their rate) are likely to be 
related to the properties of the jet such as the jet ejection rate, 
structure and kinematics of the sub-parsec scale jet. 
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Fig. 9. A sketch of the nuclear region in 3C 390.3. The draw¬ 
ing is made not to scale and shows only the approaching jet. 
The bro ad-line emission is l ikely to be generated both near 
the disk jPeterson et al.l2002l) ( BLR 1. ionized by the emission 
from a hot corona ( lFabianl2004t|Ponti et al.l2004) or the accre¬ 
tion dis k; Field & Rogers 1993) and in a rotatin g subrelativistic 
outflow dMurrav et alJ 1997HProgaetalJ2nnnl) surrounding the 
jet (BLR 2, ionized by the emission from the relativistic plasma 
in the jet). 


The results presented above imply strongly that the bulk of 
variable optical continuum and the broad-line emission asso¬ 
ciated with it are generated at a large distance from the cen¬ 
tral engine (in addition to the likely contributions from the 
hidden, “conventional” BLR close to the nucleus). Adopting 
an angle^ of 50° between the jet and the line of the sight, 
the corresponding de-projected distance between D and SI is 
0.3pc/sin50° 0.4pc. The optical continuum emission of 
3C 390.3 l eads the double-pe^ed HB line emission by ~ 30- 
100 days dShanovalova et al1l200lt Sergeev et al.ll2003i . The 
localization of the source of optical continuum with the inner¬ 
most part of the jet near S1 implies that the broad line emis¬ 
sion originates in the cone within ~ 100 light days at a dis¬ 
tance ^ 0.4 pc from the central engine. The conical shape of 
the BLR associated with the outflow arises from a large frac¬ 
tion of the continuum radiation being beamed into a cone with 
a half-opening angle^ ij/ « 30°. 

The structure of the nuclear region implied by the results 
presented above is sketched in Fig.|^ In 3C 390.3, the changes 
in the blue and red wings of Hg emiss ion line occur quasi- 
simultaneously dShanovalova et alJ200lli implying that broad¬ 
line region associated with the component S1 (outflow related 
BLR 2 in Fig. is shaping both the blue and red wings (oth¬ 
erwise one should expect significant time delay between the 
variability of the wings). The observed double-peaked struc¬ 
ture of the H/1 line can be reproduced by a rotation of the ap- 
proachingjet and a subrelativistic outflow dMurrav et al Jl 99^ 


^ On the assumption that the pattern speed and hulk speed of the jet 
of 3C 390.3 are equal, the jet inclination angle 9 ~ 50°, hulk Lorentz 
factor y ~ 2 and beaming angle <1/ x ~ 30° are estimated usin g the 
variable Doppler factor 6 = 1.16 dLahteenmaki & Valtaoiall999l) and 
the maximum apparent speed of 1.5 c observed in the compact jet. 


IProga et~ lJ200(lll . BLR 2 is evident in the broad-line emission 
at times when the jet emission dominates the optical continuum 
(see Fig. 13 }. BLR 1 may be manifested in the broad-line emis¬ 
sion around the epochs of minima in the continuum flux, when 
the jet contribution to the ionizing continuum is small. 

4. Discussion and conciusions 

The large distance of the BLR 2 from the central engine chal¬ 
lenges the existing models in which the broad-line emission is 
local ized exclusively arou nd the disk or near the central en¬ 
gine dPeterson et alfcoQ^l . It also questions the assumption of 
virialized motion in the BLR dKasni et alJl2000l) . which forms 
the foundation of the method for estimating black h ole masses 
from reverberation mapping dPeterson et alJl2nn2ll . Time de¬ 
lays and profile widths measured during periods when the jet 
emission is dominant may not necessarily reflect the Keplerian 
motion in the disk, but rather trace the rotation and outward 
motion in an outflow. This can result in large errors in esti¬ 
mates of black hole masses made from monitoring of the broad¬ 
line emission. In the case of 3C 390.3, the black hole mass 
(2.1 X 10® solar masses, Mq) estimated effectively from the 
measurements near the m aximum in the continuum light curve 
dShanovalova et alJ200ll) is significantly l arger than the values 
(3.5^ X 10^ Mn) reported in other works ( IWandel et alJll99^ 
iKasni et alJ200fll) . This difference is reconciled by considering 
the line width and the time delay between the optical contin¬ 
uum and line fluxes near the minimum of the continuum light 
curve, which yields Mbh = 3.8 x 10^ Mq. The possible exis¬ 
tence of an outflow-like region in a number of radio-loud AGN 
should be taken into account when estimates of the nuclear 
mass are made from the variability of broad emission lines. 

While the continuum emission at the base of the counterjet 
is likely to be too weak to be detected because of the relativis¬ 
tic dimming and large opacity in the disk, the optical/UV/X- 
ray line emissions from BLR 2 in the counterjet are free from 
relativistic effects and have better chances to penetrate the ab¬ 
sorbing medium towards the observer. Detection of a time lag 
between correlated variability of emission lines from the far- 
side and near-side BLRs would be the most direct observational 
evidence for double BLRs ionized by the continuum radiation 
from the bases of the jet and counterjet. If the jets are intrinsi¬ 
cally symmetric then a time delay of ~ 2 years is expected for 
such variations in 3C 390.3 (assuming the distance of ~ 1 pc 
between oppositely positioned BLRs and the jet viewing angle 
of 50°). 

The presence of the BLR2 in radio-loud AGN is capable 
of explaining some of the spectral characteristics of emission 
lines. Depending on the orientation of the jet, the approaching 
and rotating outflow material in the BLR2 will imprint promi¬ 
nent signatures on the emission lines. At small viewing angles 
of the jet the BLR2 may produce blueshifted and single-peaked 
broad emi ssion lines, while non-shifte d and double-peak emis¬ 
sion lines jEracleous & HalDerrJ2003() will be observed at large 
angles of the BLR2 to the line of si ght. Similar cha racteristics 
will have the narrow emission lines llBorosorl2003i ionized in 
the approaching subrelativistic outflow by the beamed contin¬ 
uum emission of the jet. 
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The principal results of this work and its implications are: 

1. Analysis of combined radio VLBI, optical/UV and X- 
ray data reveals significant correlations between variable 
optical continuum flux (5100A) and radio flux density 
(15GHz) of D and SI stationary components of the jet 
(at the 99% and 95% confidence levels respectively). The 
optical emission follows the radio flares with the mean 
To-opt ~ 1.4 yr and Tsi-opt ~ 0.4 yr. This finding indicates 
for the physical link between the jet and optical continuum: 
the variable optical continuum emission is located in the 
innermost part of the jet near to component SI and it is of 
non-thermal origin. This link is also supported by the cor¬ 
relation between the local maxima in the optical continuum 
light curve and the epochs at which the moving components 
of the jet pass the stationary radio feature S1. 

2. These results have important implications for the structure 
of the sub-parsec-scale nuclear region. We suggest that (i) 
the characteristics of the long-term variability of optical 
continuum emission are related to the properties of the sub¬ 
parsec scale jet (ejection rate of the radio components, its 
structure and kinematics), (ii) the bulk of optical variable 
emission originates in the jet at a distance more than or 
equal to 0.4 pc from the central engine, and (iii) the contin¬ 
uum radiation of the jet forms two BLRs (associated with 
the jet and counterjet) in the subrelativistic outflow around 
the jet. 
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APPENDIX: Properties of the VLBA images and 
model fits of 3C 390.3 

Table^compares the CLEAN component models (denoted “I”) 
and Gaussian model fits (denoted “M”) for the VLBA data 
used in the paper. The columns are: S total - total flux density 
[mJy/beam]; S peak - peak flux density [mJy/beam]; 5 mm - min¬ 
imum flux density [mJy/beam]; - goodness of the fit pa¬ 
rameter; CTjms.uv - root-mean-square between the observed and 
model visibilities [mJy]. Last row presents average ratios be¬ 
tween the respective image and model fit properties. The ra¬ 
tios are close to unity for the total and peak flux densities. The 
parameter and the visibility r.m.s. are only slightly larger 
for the Gaussian model fits, which indicates that the model fits 
represent the structure adequately. The higher maximum nega¬ 
tive flux density in the model fits (column 5 min) indicates that 
that the SNR of the Gaussian fits is on average 1.5 times lower 
than that of the VLBI images (due to increased non-Gaussian 
shapes of low-brightness regions). This reduction does not af¬ 
fect the fitted values of the component parameters (albeit it does 
increase the parameter errors) as it is related to extended emis¬ 
sion associated with the underlying flow. 

List of Objects 

3C^9R^ on page 2 

■3C 12nf on page 5 
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Table 1. Properties of the VLBA images and model fits of 3C 390.3. 



1994.67 423.9 422.9 213.3 211.1 -2.2 -3.4 1.175 1.186 253.1 253.8 

1995.96 463.1 466.2 190.8 190.8 -1.6 -2.8 0.328 0.334 119.2 120.3 

1996.37 471.6 477.6 192.4 192.6 -1.3 -3.0 0.453 0.475 83.4 85.5 

1996.82 389.5 399.1 174.2 174.2 -1.5 -3.2 1.235 1.504 42.5 45.2 

1997.19 370.0 377.5 179.1 180.4 -0.7 -2.2 1.533 1.800 39.9 43.0 

1997.66 418.3 425.2 267.4 267.2 -1.8 -2.2 0.575 0.580 152.1 152.5 

1998.21 392.4 389.6 191.2 191.4 -1.1 -1.3 0.429 0.430 94.6 94.7 

1999.85 313.2 315.8 216.9 218.6 -1.1 -1.3 1.103 1.149 31.0 31.4 

2001.17 274.6 274.3 203.6 204.2 -1.1 -1.2 1.211 1.216 30.7 30.8 

2002.13 308.8 303.6 208.4 210.3 -0.9 -1.2 1.020 1.146 26.5 27.8 

< I/M > 0.994 0.998 0.664 0.943 0.977 














